Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,382,000 C→T 8.5% C18C (TGC→TGT ompG → outer membrane porin G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,382,0000CT 8.5% 290.1 / 14.0 94C18C (TGC→TGTompGouter membrane porin G
Reads supporting (aligned to +/- strand):  ref base C (44/42);  new base T (4/4);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

TTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTG  >  NC_000913/1381866‑1382145
                                                                                                                                      |                                                                                                                                                 
tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagga                                                                                                                                     >  4:667516/1‑149 (MQ=255)
  tCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                   <  3:667516/149‑1 (MQ=255)
    tttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagg                                                                                                                                 <  3:181853/149‑1 (MQ=255)
        aTGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagg                                                                                                                                 <  4:1188947/145‑1 (MQ=255)
        aTGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagg                                                                                                                                 >  3:1188947/1‑145 (MQ=255)
               aaCGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAGCAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGcc                                                                                                                                          <  1:1248023/129‑1 (MQ=255)
               aaCGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGcc                                                                                                                                          >  2:1248023/1‑129 (MQ=255)
                     aaTAGCAATTCGCTAAATACGGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaaggaa                                                                                                                               >  2:1357629/1‑134 (MQ=255)
                     aaTAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaaggaa                                                                                                                               <  1:1357629/134‑1 (MQ=255)
                             ttCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGcca                                                                                                                                         >  2:1415542/1‑115 (MQ=255)
                             ttCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGcca                                                                                                                                         <  1:1415542/116‑2 (MQ=255)
                             ttCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcg                                                                                                        >  4:176671/1‑149 (MQ=255)
                                gCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGa                                                                                                     <  3:176671/149‑1 (MQ=255)
                                       aCAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGa                                                                                                                            >  1:1417636/1‑119 (MQ=255)
                                       aCAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGa                                                                                                                            <  2:1417636/119‑1 (MQ=255)
                                         tgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                        >  1:437081/2‑105 (MQ=255)
                                         tgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                        <  2:437081/104‑1 (MQ=255)
                                         aGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                            <  1:108386/149‑1 (MQ=255)
                                           ggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACTGGCCGaggaaa                                                                                                                                   <  3:860848/108‑1 (MQ=255)
                                           ggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                   >  4:860848/1‑108 (MQ=255)
                                           ggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                   >  3:293591/1‑108 (MQ=255)
                                           ggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                   <  4:293591/108‑1 (MQ=255)
                                           ggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGag                                                                                                                                       <  1:1203727/104‑1 (MQ=255)
                                           ggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGag                                                                                                                                       >  2:1203727/1‑104 (MQ=255)
                                                   cATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAAc                                                                                                                              >  2:152356/1‑105 (MQ=255)
                                                   cATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAAc                                                                                                                              <  1:152356/105‑1 (MQ=255)
                                                    aTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAg                                                                                                                                             <  3:661372/89‑1 (MQ=255)
                                                    aTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAg                                                                                                                                             >  4:661372/1‑89 (MQ=255)
                                                         cccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                   >  4:126062/1‑94 (MQ=255)
                                                         cccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                   >  1:1228998/1‑94 (MQ=255)
                                                         cccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                   <  2:1228998/94‑1 (MQ=255)
                                                         cccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                   <  3:126062/94‑1 (MQ=255)
                                                         cccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAata                                                                                                              >  2:449416/1‑115 (MQ=255)
                                                         cccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAata                                                                                                              <  1:449416/115‑1 (MQ=255)
                                                          ccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                            <  4:568464/132‑1 (MQ=255)
                                                          ccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                            >  3:568464/1‑132 (MQ=255)
                                                                   taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTg                                                                                                                         <  3:1309705/94‑1 (MQ=255)
                                                                   taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTg                                                                                                                         >  4:1309705/1‑94 (MQ=255)
                                                                   taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGAt                                                                                                    >  4:247033/1‑115 (MQ=255)
                                                                   taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGAt                                                                                                    <  3:1329435/115‑1 (MQ=255)
                                                                   taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGAt                                                                                                    <  3:247033/115‑1 (MQ=255)
                                                                   taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGAt                                                                                                    >  4:1329435/1‑115 (MQ=255)
                                                                   taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGAt                                                                                                    >  2:997429/1‑115 (MQ=255)
                                                                   taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGAt                                                                                                    <  1:997429/115‑1 (MQ=255)
                                                                    aCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAatat                                                                                                             <  1:1366848/105‑1 (MQ=255)
                                                                    aCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAatat                                                                                                             >  2:1366848/1‑105 (MQ=255)
                                                                    aCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGtt                                                                          <  1:1139553/140‑1 (MQ=255)
                                                                    aCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGtt                                                                          >  2:1139553/1‑140 (MQ=255)
                                                                    aCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAg                                                                 >  3:1329780/1‑149 (MQ=255)
                                                                      aaGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGag                                                                                                                                       >  4:80929/1‑77 (MQ=255)
                                                                      aaGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGag                                                                                                                                       <  3:80929/77‑1 (MQ=255)
                                                                        ggAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcgc                                                                                                       <  4:26952/107‑1 (MQ=255)
                                                                        ggAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcgc                                                                                                       >  3:26952/1‑107 (MQ=255)
                                                                         gAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAg                                                                                                                                             <  2:346365/68‑1 (MQ=255)
                                                                         gAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAg                                                                                                                                             >  1:346365/1‑68 (MQ=255)
                                                                         gAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGt                                                                                  >  4:593407/1‑127 (MQ=255)
                                                                         gAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGt                                                                                  <  3:593407/127‑1 (MQ=255)
                                                                            cAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggat                                                         <  4:1329780/149‑1 (MQ=255)
                                                                                   gAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTgg                                                  >  4:1335132/1‑149 (MQ=255)
                                                                                    aaaaaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGAc                                                                                                                           >  1:37864/1‑75 (MQ=255)
                                                                                    aaaaaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGAc                                                                                                                           <  2:37864/75‑1 (MQ=255)
                                                                                      aaaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCgg                                               >  3:1432130/1‑149 (MQ=255)
                                                                                          ttattaCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGcc                                           >  1:693654/1‑149 (MQ=255)
                                                                                               aCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACt                                                                                                                    <  4:108237/71‑1 (MQ=255)
                                                                                               aCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACt                                                                                                                    >  3:108237/1‑71 (MQ=255)
                                                                                                cccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGt                                     >  1:376234/1‑149 (MQ=255)
                                                                                                 ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                            >  3:286229/1‑93 (MQ=255)
                                                                                                 ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                            <  4:286229/93‑1 (MQ=255)
                                                                                                 ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGtt                                                                          >  3:675751/1‑111 (MQ=255)
                                                                                                 ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGtt                                                                          <  4:675751/111‑1 (MQ=255)
                                                                                                 ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGaa                                                                  >  2:1346811/1‑119 (MQ=255)
                                                                                                 ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGaa                                                                  <  1:1346811/119‑1 (MQ=255)
                                                                                                  cTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCt                                   <  2:693654/149‑1 (MQ=255)
                                                                                                             tGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATg                                                                       >  1:232504/1‑102 (MQ=255)
                                                                                                             tGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATg                                                                       <  2:232504/102‑1 (MQ=255)
                                                                                                               gTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTaa                            >  1:893394/1‑143 (MQ=255)
                                                                                                               gTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTaa                            <  2:893394/143‑1 (MQ=255)
                                                                                                                    gtgtgCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTCATATCGGCGCGATGTACGAAATAGAAAACGTCGa                                                                               <  3:1293036/87‑1 (MQ=255)
                                                                                                                    gtgtgCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGa                                                                               >  4:1293036/1‑87 (MQ=255)
                                                                                                                         cGGGAATTGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggat                                                         <  2:1424904/104‑1 (MQ=255)
                                                                                                                         cGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATg                                                                       >  4:1031327/1‑90 (MQ=255)
                                                                                                                         cGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATg                                                                       <  3:1031327/90‑1 (MQ=255)
                                                                                                                         cGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggat                                                         >  1:1424904/1‑104 (MQ=255)
                                                                                                                          gggAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGt           >  2:1328827/1‑149 (MQ=255)
                                                                                                                           ggAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTg          >  1:1216926/1‑149 (MQ=255)
                                                                                                                            gAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTgg         >  2:1388105/1‑149 (MQ=255)
                                                                                                                            gAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTgg         <  3:1335132/149‑1 (MQ=255)
                                                                                                                              tcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGc                                                                     <  4:163552/85‑1 (MQ=255)
                                                                                                                              tcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGc                                                                     >  3:163552/3‑87 (MQ=255)
                                                                                                                              tcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGa                                                                >  1:719953/3‑92 (MQ=255)
                                                                                                                              tcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGa                                                                <  2:719953/90‑1 (MQ=255)
                                                                                                                                 gCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAAt    <  2:376234/149‑1 (MQ=255)
                                                                                                                                   cTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTg  <  2:1216926/149‑1 (MQ=255)
                                                                                                                                   cTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTg  >  2:600868/1‑149 (MQ=255)
                                                                                                                                      |                                                                                                                                                 
TTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTG  >  NC_000913/1381866‑1382145

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: