Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,993,569 C→T 100% P806L (CCG→CTG)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,993,5690CT97.8% 770.4 / 2.0 208P806L (CCG→CTG) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (107/96);  total (109/99)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 210 total aligned reads displayed.

ATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAA  >  NC_000913/3993421‑3993713
                                                                                                                                                    |                                                                                                                                                
atCTACATTATCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGTCAGCAAAGAGGAGCTGATACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGcc                                                                                                                                                  >  3:940486/1‑149 (MQ=255)
   tACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCa                                                                                                                                               >  2:924729/1‑149 (MQ=255)
    aCATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAg                                                                                                                                              >  3:311431/1‑149 (MQ=255)
     cATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGt                                                                                                                                             <  3:116427/149‑1 (MQ=255)
        tctcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCt                                                                                                                                          <  4:1112849/149‑1 (MQ=255)
          tcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTAt                                                                                                                                        >  4:782686/1‑149 (MQ=255)
          tcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTAt                                                                                                                                        <  1:924729/149‑1 (MQ=255)
          tcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTAt                                                                                                                                        >  1:870958/1‑149 (MQ=255)
              cgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGa                                                                                                                                    <  3:136381/149‑1 (MQ=255)
               gaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGAt                                                                                                                                   >  3:1296909/1‑149 (MQ=255)
                  aGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATtgt                                                                                                                                >  2:504732/1‑149 (MQ=255)
                     aaCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATc                                                                                                                                       <  1:117924/139‑1 (MQ=255)
                     aaCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATc                                                                                                                                       >  2:117924/1‑139 (MQ=255)
                     aaCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGAtt                                                                                                                                  <  1:83789/144‑1 (MQ=255)
                      aCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAg                                                                                                                            >  1:944611/1‑149 (MQ=255)
                       ccGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATtgt                                                                                                                                <  1:1403674/144‑1 (MQ=255)
                          ggTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTg                                                                                                                        >  1:433933/1‑149 (MQ=255)
                           gTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTa                                                                                                                                         >  4:798231/1‑131 (MQ=255)
                            ttGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGAt                                                                                                                      <  1:1461975/149‑1 (MQ=255)
                            ttGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGAt                                                                                                                      <  4:311431/149‑1 (MQ=255)
                                ggTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTAt                                                                                                                                        <  2:1274921/127‑1 (MQ=255)
                                       caccacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACa                                                                                                           >  2:605227/1‑149 (MQ=255)
                                        accacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAg                                                                                                          >  2:635044/1‑149 (MQ=255)
                                         ccacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAgg                                                                                                         >  2:1016935/1‑149 (MQ=255)
                                          cacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGt                                                                                                        >  1:396685/1‑149 (MQ=255)
                                             tGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAt                                                                                                     <  4:1296909/149‑1 (MQ=255)
                                             tGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAt                                                                                                     >  2:1134733/1‑149 (MQ=255)
                                               cGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGaa                                                                                                                             <  3:38032/123‑1 (MQ=255)
                                               cGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttc                                                                                                   <  1:504732/149‑1 (MQ=255)
                                               cGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttc                                                                                                   <  2:396685/149‑1 (MQ=255)
                                               cGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttc                                                                                                   >  3:787350/1‑149 (MQ=255)
                                                       gcaAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGc                                                                                           <  1:605227/149‑1 (MQ=255)
                                                         aaaGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGGTGGTCGTGAACAGGTGAtt                                                                                                    >  1:28631/1‑138 (MQ=255)
                                                          aaGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATc                                                                                                                                       <  3:476887/102‑1 (MQ=255)
                                                             aggagCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACaaaa                                                                                     <  2:433933/149‑1 (MQ=255)
                                                                tgCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGt                                                                                                                >  3:434257/2‑119 (MQ=255)
                                                                     gTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTAt                                                                                >  4:7923/1‑146 (MQ=255)
                                                                      tACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGt                                                                                                                <  4:391633/113‑1 (MQ=255)
                                                                         gTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTAt                                                                                >  1:970199/1‑142 (MQ=255)
                                                                          tGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGAtt                                                                                                                                  >  1:1311777/1‑91 (MQ=255)
                                                                           gACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAgg                                                                                                         <  1:1130465/115‑1 (MQ=255)
                                                                               tCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACaaaa                                                                                     <  1:735174/131‑1 (MQ=255)
                                                                               tCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACaaaa                                                                                     >  2:735174/1‑131 (MQ=255)
                                                                                   tCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATtgtg                                                                                                                               <  4:662298/85‑1 (MQ=255)
                                                                                   tCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATtgtg                                                                                                                               >  3:662298/1‑85 (MQ=255)
                                                                                       ttCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACa                                                                                                           <  4:364800/101‑1 (MQ=255)
                                                                                          tACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACa                                                                                                           >  2:843325/1‑98 (MQ=255)
                                                                                            cTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGAt                                                                                                                      <  2:1083444/85‑1 (MQ=255)
                                                                                            cTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTg                                                               <  4:714243/140‑1 (MQ=255)
                                                                                             tCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTAt                                                                                >  2:820126/1‑122 (MQ=255)
                                                                                             tCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGa                                                     >  3:713779/1‑149 (MQ=255)
                                                                                               gTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTaa                                                                          <  2:64389/126‑1 (MQ=255)
                                                                                               gTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATc                                                   >  1:1256989/1‑149 (MQ=255)
                                                                                               gTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATc                                                   >  3:350118/1‑149 (MQ=255)
                                                                                                tCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACaaaa                                                                                     >  3:172096/1‑114 (MQ=255)
                                                                                                  aTCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGa                                                                                                                       <  2:349534/78‑1 (MQ=255)
                                                                                                  aTCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGc                                                              <  2:969513/135‑1 (MQ=255)
                                                                                                  aTCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACg                                                >  3:794686/1‑149 (MQ=255)
                                                                                                   tCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGa                                               <  3:852064/149‑1 (MQ=255)
                                                                                                   tCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGa                                               >  1:1078400/1‑149 (MQ=255)
                                                                                                     gCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTAt                                                                                <  3:360255/114‑1 (MQ=255)
                                                                                                     gCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATa                                             <  1:1134733/149‑1 (MQ=255)
                                                                                                          aCCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAg                                                       <  3:593819/134‑1 (MQ=255)
                                                                                                          aCCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAg                                                       >  4:593819/1‑134 (MQ=255)
                                                                                                           ccGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAcgccgc                                       >  2:1414387/1‑149 (MQ=255)
                                                                                                           ccGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAc                                            <  4:1183592/144‑1 (MQ=255)
                                                                                                           ccGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAc                                            >  3:1183592/1‑144 (MQ=255)
                                                                                                             gTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGTTTGTGaa                                                                                                                             >  4:457746/1‑61 (MQ=255)
                                                                                                              ttttACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATgccg                                                                   >  3:753263/1‑118 (MQ=255)
                                                                                                              ttttACCTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTAt                                    >  3:816364/1‑149 (MQ=255)
                                                                                                              ttttACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTAt                                    <  4:1413368/149‑1 (MQ=255)
                                                                                                                    cTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTAt                                                                                <  3:526037/99‑1 (MQ=255)
                                                                                                                    cTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTAt                                                                                >  4:526037/1‑99 (MQ=255)
                                                                                                                    cTACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGc                              >  4:74890/1‑149 (MQ=255)
                                                                                                                     tACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGaa                                                                                                                             >  1:643171/1‑53 (MQ=255)
                                                                                                                     tACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGt                                                                                                        >  4:1407745/1‑74 (MQ=255)
                                                                                                                     tACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGt                                                                                                        <  3:1407745/74‑1 (MQ=255)
                                                                                                                     tACGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCa                             >  2:169788/1‑149 (MQ=255)
                                                                                                                      aCGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTa                                  <  1:680571/143‑1 (MQ=255)
                                                                                                                      aCGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTa                                  >  2:680571/1‑143 (MQ=255)
                                                                                                                      aCGGCTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCaa                            >  3:208978/1‑149 (MQ=255)
                                                                                                                          cTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCAc                                                 <  4:625475/124‑1 (MQ=255)
                                                                                                                          cTCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCAc                                                 >  3:625475/1‑124 (MQ=255)
                                                                                                                           tCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTAt                                                                                <  1:610948/92‑1 (MQ=255)
                                                                                                                           tCAAGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATAttt                       >  4:1020430/1‑149 (MQ=255)
                                                                                                                             aaGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATAttttt                     >  4:1253078/1‑149 (MQ=255)
                                                                                                                             aaGCTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATAttttt                     >  3:1410750/1‑149 (MQ=255)
                                                                                                                                cTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGt                  <  3:744528/149‑1 (MQ=255)
                                                                                                                                cTTCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGt                  <  3:1020430/149‑1 (MQ=255)
                                                                                                                                 ttCATCAACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtg                 >  2:240155/1‑149 (MQ=255)
                                                                                                                                  tcatcaACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtga                <  3:1253078/149‑1 (MQ=255)
                                                                                                                                     tcaACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgatga             >  1:72064/1‑149 (MQ=255)
                                                                                                                                     tcaACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgatga             <  4:208978/149‑1 (MQ=255)
                                                                                                                                     tcaACTTCAACCTGCTGCAGTTCTATCAGATCGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgatga             >  1:650303/1‑149 (MQ=255)
                                                                                                                                      caACTTCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGaa            <  4:350118/149‑1 (MQ=255)
                                                                                                                                          ttCAACCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgat               >  1:100557/1‑142 (MQ=255)
                                                                                                                                             aaCCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATAtttt                      >  2:1316625/1‑132 (MQ=255)
                                                                                                                                              aCCTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAcgc                                          >  2:33065/1‑111 (MQ=255)
                                                                                                                                               ccTGCTGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGa   >  4:872152/1‑149 (MQ=255)
                                                                                                                                               ccTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGa   >  2:72520/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAA  >  NC_000913/3993421‑3993713

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