Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 236,388 C→G 5.5% L149L (CTG→CTC argH ← argininosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947236,3880CG 5.5% 285.4 / 5.0 90L149L (CTG→CTCargHargininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base C (46/39);  new base G (3/2);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GCGCTTGCGGCAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCGGTCGGTCAGGCG  >  NC_002947/236241‑236537
                                                                                                                                                   |                                                                                                                                                      
gcgcTTGCGGCAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCcag                                                                                                                                                      >  1:397935/1‑150 (MQ=255)
gcgcTTGCGGCAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCcag                                                                                                                                                      >  2:2409891/1‑150 (MQ=255)
gcgcTTGCGGCAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCcag                                                                                                                                                      >  1:418828/1‑150 (MQ=255)
       cGGCAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGccc                                                                                                                                               >  1:1211964/1‑150 (MQ=255)
          cAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCctgc                                                                                                                                            <  2:2163439/150‑1 (MQ=255)
          cAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCctgc                                                                                                                                            <  2:286662/150‑1 (MQ=255)
          cAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCCGCCTGCTGGAGCAGGcagagc                                                                                                                                            <  1:1017177/150‑6 (MQ=255)
            gTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGTTGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCctgctg                                                                                                                                          >  1:1071455/1‑150 (MQ=255)
                aCCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAgg                                                                                                                                      >  1:935230/1‑150 (MQ=255)
                 ccAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGc                                                                                                                                     >  1:882337/1‑150 (MQ=255)
                  cAGGCGCTCATAGTCGCAACAGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGTTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCg                                                                                                                                    >  2:217051/1‑150 (MQ=255)
                   aGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCTGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCgg                                                                                                                                   <  1:1209903/150‑1 (MQ=255)
                         tCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTTCCTGCTCCAGCAGGCCCTGCTGCAGGGGGGTGAtt                                                                                                                             >  2:1142782/1‑150 (MQ=255)
                         tCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCgg                                                                                                                                   >  1:1585356/1‑144 (MQ=255)
                         tCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCgg                                                                                                                                   <  2:1585356/144‑1 (MQ=255)
                         tCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGAtt                                                                                                                             <  2:418828/150‑1 (MQ=255)
                         tCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTACAGCAGGCCCTGCTGCAGGCGGGTGAtt                                                                                                                             >  1:2177331/1‑150 (MQ=255)
                           aTAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTc                                                                                                                           >  1:1501486/1‑150 (MQ=255)
                             ggTCGCGGCTGAGCATTTCGACCCACGCCAGCAGGTGGTGGCAAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGc                                                                                                                         <  2:1656492/149‑1 (MQ=255)
                                       gAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGAAGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAg                                                                                                              <  2:342004/150‑1 (MQ=255)
                                         gCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGTTGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGt                                                                                                            <  2:1071455/150‑1 (MQ=255)
                                         acATTGCTAACCACGCTCGTAACTGGTGGCCCACGGTCCCCGGCTGCGCGGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCATGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGt                                                                                                            <  2:397935/149‑1 (MQ=255)
                                          cATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTc                                                                                                           >  1:1493766/1‑150 (MQ=255)
                                                    cACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGc                                                                                                                         <  2:1391690/127‑1 (MQ=255)
                                                    cACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGc                                                                                                                         >  1:1391690/1‑127 (MQ=255)
                                                    cACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTtcgtcg                                                                                                 >  2:969354/1‑150 (MQ=255)
                                                     aCGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTcgc                                                                                                >  2:539343/1‑150 (MQ=255)
                                                     aCGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTcgc                                                                                                <  2:1211964/150‑1 (MQ=255)
                                                       gCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCgcgc                                                                                              >  1:76383/1‑150 (MQ=255)
                                                               gtggtgGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGccc                                                                                                                      <  2:741936/118‑1 (MQ=255)
                                                               gtggtgGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGccc                                                                                                                      >  1:741936/1‑118 (MQ=255)
                                                                tggtggCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAgg                                                                                     >  1:1930814/1‑150 (MQ=255)
                                                                  gtggCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCg                                                                                   <  1:2409891/150‑1 (MQ=255)
                                                                  gtggCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCg                                                                                   <  1:1840734/150‑1 (MQ=255)
                                                                    ggCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGa                                                                                 <  1:969354/150‑1 (MQ=255)
                                                                       ttcaGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGt                                                                              >  1:85022/4‑150 (MQ=255)
                                                                       cAAAGGTGACCGGCTGCGCCGTCTGGAGGTGGGTCAAACCAGGCATGATGAGGTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGCTCGTTGCCGTCGCGCGGCGTGAGGTGGATgg                                                                              >  2:338996/1‑149 (MQ=255)
                                                                       cAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGc                                                                                                                         >  2:281671/1‑108 (MQ=255)
                                                                       cAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGc                                                                                                                         <  1:281671/108‑1 (MQ=255)
                                                                       cAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCGTGACGACGGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGt                                                                              >  1:1330038/1‑150 (MQ=255)
                                                                       cAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCgg                                                                                  <  1:1803528/146‑1 (MQ=255)
                                                                       cAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCgg                                                                                  >  2:1803528/1‑146 (MQ=255)
                                                                       cAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGt                                                                              >  1:301568/1‑150 (MQ=255)
                                                                        aaaGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTc                                                                             >  1:1136896/1‑150 (MQ=255)
                                                                        aaaGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGCCTTCACGCTCTGCCTGCACCAGTAGGCCCTGCTCCCGCCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCGGCCAACGGCCGATGTc                                                                             >  2:1068442/1‑150 (MQ=255)
                                                                           ggTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGt                                                                          <  1:948249/150‑1 (MQ=255)
                                                                           ggTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGt                                                                          >  2:143237/1‑150 (MQ=255)
                                                                           ggTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGt                                                                          <  1:104990/150‑1 (MQ=255)
                                                                           ggTGACCGGCTGCGACGTCTGCAGGTGGGTGGACCCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTTCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGt                                                                          <  2:1330038/150‑1 (MQ=255)
                                                                                  ggCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCcagc                                                                                                                                                     >  2:1014124/1‑69 (MQ=255)
                                                                                  ggCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCcagc                                                                                                                                                     <  1:1014124/69‑1 (MQ=255)
                                                                                       cgcCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTc                                                              >  2:463797/1‑150 (MQ=255)
                                                                                        gcCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGt                                                                          <  2:918579/137‑1 (MQ=255)
                                                                                        gcCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGt                                                                          >  1:918579/1‑137 (MQ=255)
                                                                                         ccGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGAtt                                                                                                                             <  2:2123256/86‑1 (MQ=255)
                                                                                         ccGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGAtt                                                                                                                             >  1:2123256/1‑86 (MQ=255)
                                                                                         ccGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGt                                                            <  2:2177331/150‑1 (MQ=255)
                                                                                          cggcggCAGGTGGCTGAAACCAGGCATGAGGGTTTGGGCTTCACCCTCTGCCTGCTCCAGCAGGCCCTGGTGCAGGCGGGTGATTTCGCCCAGGAGCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGtt                                                           <  2:973946/145‑1 (MQ=255)
                                                                                               gCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGc                                                      <  2:1930814/150‑1 (MQ=255)
                                                                                                   gtgggtgAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGcc                                                                                                                       >  1:1442772/1‑81 (MQ=255)
                                                                                                   gtgggtgAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGcc                                                                                                                       <  2:1442772/81‑1 (MQ=255)
                                                                                                   gtgggtgAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACg                                                  >  2:105909/1‑150 (MQ=255)
                                                                                                    tgggtgAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCtgctgc                                                                                                                                         >  2:261444/1‑63 (MQ=255)
                                                                                                    tgggtgAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCtgctgc                                                                                                                                         <  1:261444/63‑1 (MQ=255)
                                                                                                       gtgAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAATGCGGATGTCGGTGGCCACCTGGt                                                               <  1:3400/133‑1 (MQ=255)
                                                                                                       gtgAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAATGCGGATGTCGGTGGCCACCTGGt                                                               >  2:3400/1‑133 (MQ=255)
                                                                                                       gtgAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCa                                              <  1:478948/150‑1 (MQ=255)
                                                                                                                 gCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGc                                                         >  1:838945/1‑129 (MQ=255)
                                                                                                                 gCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGc                                                         <  2:838945/129‑1 (MQ=255)
                                                                                                                     gatgGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGcttctt                                >  2:984282/1‑150 (MQ=255)
                                                                                                                     gatgGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGcttctt                                <  2:133856/150‑1 (MQ=255)
                                                                                                                      atgGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGcttctt                                >  2:778895/1‑149 (MQ=255)
                                                                                                                      atgGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGcttctt                                <  1:778895/149‑1 (MQ=255)
                                                                                                                       tgGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGc                              <  2:76383/150‑1 (MQ=255)
                                                                                                                         gTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCg                            >  2:2263074/1‑150 (MQ=255)
                                                                                                                         gTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCg                            >  1:554125/1‑150 (MQ=255)
                                                                                                                          tttCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAGAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCgg                           >  1:603556/1‑150 (MQ=255)
                                                                                                                                 ttCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCa                                                                     <  1:2352990/101‑1 (MQ=255)
                                                                                                                                 ttCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCa                                                                     >  2:2352990/1‑101 (MQ=255)
                                                                                                                                    aCGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGAt                 <  1:2263074/150‑1 (MQ=255)
                                                                                                                                      gCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGc               >  2:2258452/1‑150 (MQ=255)
                                                                                                                                        tctGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGcgc             <  2:1649338/150‑2 (MQ=255)
                                                                                                                                         gggcctgcTGGAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGc                                                                                    >  1:1064066/3‑78 (MQ=255)
                                                                                                                                         gggcctgcTGGAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGc                                                                                    <  2:1064066/76‑1 (MQ=255)
                                                                                                                                         gggcctgcTGGAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGcggt            >  2:1248480/3‑150 (MQ=255)
                                                                                                                                          ggcctgcTGGAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGcggtc           >  2:2336371/2‑150 (MQ=255)
                                                                                                                                           gcctgcTCCAGCAATCCCTTCTGAAGGCGCGTGATTTCGCCCAGGATCAGGTCGATTTCGTCTCCCATCCAAAGGCGTATGTCGGTCGCCACCTTTTCGTTGCAGCTACGCTCAGTATTCAGCTTCTTGCCCTTTATGCCGATGcggtcg          <  2:178210/150‑1 (MQ=255)
                                                                                                                                            cctgcTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCAcc                                                                   <  2:1315135/92‑1 (MQ=255)
                                                                                                                                            cctgcTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCAcc                                                                   >  1:1315135/1‑92 (MQ=255)
                                                                                                                                            cctgcTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGcggtcgg         >  2:2326265/1‑150 (MQ=255)
                                                                                                                                            cctgcTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCAGGTGATGCCGATGcggtcgg         >  1:1291689/1‑150 (MQ=255)
                                                                                                                                            ccgccTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGcggtcgg         >  1:1493895/5‑150 (MQ=255)
                                                                                                                                              tgcTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCggtcggtc       >  2:2050015/1‑150 (MQ=255)
                                                                                                                                                 tCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGtt                                                           >  2:1992518/1‑95 (MQ=255)
                                                                                                                                                 tCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGtt                                                           <  1:1992518/95‑1 (MQ=255)
                                                                                                                                                   cagcagGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCa                                       <  2:347863/113‑1 (MQ=255)
                                                                                                                                                   cagcagGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCa                                       >  1:347863/1‑113 (MQ=255)
                                                                                                                                                   cagcagGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCGGTCGGTCAGGCg  >  1:707797/1‑150 (MQ=255)
                                                                                                                                                   cagcagGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCGGTCGGTCAGGCg  >  2:2095672/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GCGCTTGCGGCAGTCGACCAGGCGCTCATAGTCGCGGCTGAGCATTTCGAACCACGCCAGCAGGTGGTGGCCAAAGGTGACCGGCTGCGCCGTCTGCAGGTGGGTGAAACCAGGCATGATGGTTTCGGCTTCACGCTCTGCCTGCTCCAGCAGGCCCTGCTGCAGGCGGGTGATTTCGCCCAGGATCAGGTCGATTTCGTCGCGCAGCCAAAGGCGGATGTCGGTGGCCACCTGGTCGTTGCGGCTACGCCCAGTATGCAGCTTCTTGCCGGTGATGCCGATGCGGTCGGTCAGGCG  >  NC_002947/236241‑236537

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