Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,081,426 C→A 6.7% V145V (GTG→GTT fdoI ← formate dehydrogenase‑O, cytochrome b556 subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,081,4260CA 6.7% 312.0 / 6.3 89V145V (GTG→GTTfdoIformate dehydrogenase‑O, cytochrome b556 subunit
Reads supporting (aligned to +/- strand):  ref base C (58/25);  new base A (3/3);  total (61/28)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

GCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCGGTGTCACCTACTT  >  NC_000913/4081289‑4081562
                                                                                                                                         |                                                                                                                                        
gCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCAc                                                                                                                                        <  1:1017806/140‑1 (MQ=255)
 cccAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCAcc                                                                                                                                       <  1:1559140/140‑1 (MQ=255)
         cTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTga                                                                                                                               >  2:690576/1‑140 (MQ=255)
           ggTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTgaga                                                                                                                             >  1:1298935/1‑140 (MQ=255)
           ggTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTgaga                                                                                                                             >  2:496001/1‑140 (MQ=255)
            gTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCCCCGGGATTGAGaa                                                                                                                            >  2:1271765/1‑140 (MQ=255)
                cccATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCa                                                                                                                        >  1:1167172/1‑140 (MQ=255)
                    tCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGgcg                                                                                                                    >  1:679979/1‑140 (MQ=255)
                    tCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGgcg                                                                                                                    >  1:679931/1‑140 (MQ=255)
                    tCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGgcg                                                                                                                    >  1:679910/1‑140 (MQ=255)
                     ccTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGgcgc                                                                                                                   >  2:856974/1‑140 (MQ=255)
                      cTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCa                                                                                                                  >  1:874286/1‑140 (MQ=255)
                      cTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCa                                                                                                                  >  1:412080/1‑140 (MQ=255)
                      cTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGAAAACCACAATTAACGCCACTGCGGCAAAGGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCa                                                                                                                  >  1:223237/1‑140 (MQ=255)
                          caccaTGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAAt                                                                                                              >  2:532339/1‑140 (MQ=255)
                             caTGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAg                                                                                                           <  1:902620/140‑1 (MQ=255)
                               tGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGa                                                                                                         >  2:956331/1‑140 (MQ=255)
                               tGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGa                                                                                                         >  1:293672/1‑140 (MQ=255)
                                           tGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATc                                                                                             >  2:626561/1‑140 (MQ=255)
                                           tGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATAACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATc                                                                                             >  1:636979/1‑140 (MQ=255)
                                           tGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATAACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATc                                                                                             >  1:636966/1‑140 (MQ=255)
                                           tGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATAACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATc                                                                                             >  1:636943/1‑140 (MQ=255)
                                            gCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCa                                                                                            >  2:1599326/1‑140 (MQ=255)
                                               tttCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCAcgc                                                                                         <  2:412080/140‑1 (MQ=255)
                                               tttCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCAcgc                                                                                         <  1:1579659/140‑1 (MQ=255)
                                                 tCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCAcgccg                                                                                       >  2:401620/1‑140 (MQ=255)
                                                 tCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCAcgccg                                                                                       <  1:1480920/140‑1 (MQ=255)
                                                 tCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCAcgccg                                                                                       <  1:1186871/140‑1 (MQ=255)
                                                  cACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCAcgccgc                                                                                      <  2:936496/140‑1 (MQ=255)
                                                            gcggcgTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACcagcag                                                                            <  2:1349563/140‑1 (MQ=255)
                                                             cggcgTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACcagcagc                                                                           >  2:1562484/1‑140 (MQ=255)
                                                               gcgTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACcagcagcag                                                                         >  2:1417683/1‑140 (MQ=255)
                                                               gcgTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACcagcagcag                                                                         >  2:4883/1‑140 (MQ=255)
                                                               gcgTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACcagcagcag                                                                         >  2:1417936/1‑140 (MQ=255)
                                                               gcgTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACcagcagcag                                                                         >  1:1357603/1‑140 (MQ=255)
                                                                cgTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGa                                                                        <  2:1580717/140‑1 (MQ=255)
                                                                          cACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGaaaa                                                              >  2:685290/1‑140 (MQ=255)
                                                                          cACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGaaaa                                                              >  1:676731/1‑140 (MQ=255)
                                                                          cACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGaaaa                                                              >  1:389120/1‑140 (MQ=255)
                                                                          cACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGaaaa                                                              >  1:819284/1‑140 (MQ=255)
                                                                          cACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGaaaa                                                              >  2:173549/1‑140 (MQ=255)
                                                                           aCCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAaat                                                             >  1:616538/1‑140 (MQ=255)
                                                                           aCCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAaat                                                             >  1:340276/1‑140 (MQ=255)
                                                                           aCCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAaat                                                             >  1:616594/1‑140 (MQ=255)
                                                                               tGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATc                                                         >  1:1632639/1‑140 (MQ=255)
                                                                               tGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATc                                                         >  2:316488/1‑140 (MQ=255)
                                                                               tGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATc                                                         >  2:48456/1‑140 (MQ=255)
                                                                               tGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATc                                                         >  2:1504011/1‑140 (MQ=255)
                                                                               tGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATc                                                         >  2:1147956/1‑140 (MQ=255)
                                                                               tGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCACAACCAGGAAAATAATc                                                         >  2:48498/1‑140 (MQ=255)
                                                                               tGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCACGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATc                                                         >  2:745325/1‑140 (MQ=255)
                                                                               tGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATAAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATc                                                         >  1:1195019/1‑140 (MQ=255)
                                                                                     ccACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCgccgcc                                                   <  2:1098540/140‑1 (MQ=255)
                                                                                         aaTTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGa                                               >  2:530592/1‑140 (MQ=255)
                                                                                         aaTTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGa                                               <  2:107487/140‑1 (MQ=255)
                                                                                         aaTTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGa                                               >  1:728706/1‑140 (MQ=255)
                                                                                         aaTTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGa                                               >  1:728712/1‑140 (MQ=255)
                                                                                         aaTTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAAGAATCGCCGCCCAGa                                               >  1:1685858/1‑140 (MQ=255)
                                                                                          aTTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGaa                                              >  1:688944/1‑140 (MQ=255)
                                                                                          aTTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGaa                                              >  1:122087/1‑140 (MQ=255)
                                                                                           ttAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGaaa                                             >  2:1237638/1‑140 (MQ=255)
                                                                                            tAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGaaaa                                            <  1:117440/140‑1 (MQ=255)
                                                                                            tAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGaaaa                                            <  1:502348/140‑1 (MQ=255)
                                                                                                   aCTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCAttt                                     >  1:466986/1‑140 (MQ=255)
                                                                                                   aCTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAAAGCAttt                                     >  1:1384534/1‑140 (MQ=255)
                                                                                                     tGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCt                                   >  2:590751/1‑140 (MQ=255)
                                                                                                        ggCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGAc                                <  1:345225/140‑1 (MQ=255)
                                                                                                           aaaTGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGa                             <  2:1216581/140‑1 (MQ=255)
                                                                                                             atgaatgCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAg                           <  2:1534267/140‑1 (MQ=255)
                                                                                                                  tgCAGCATTAACGCGAATCGGATAACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTtata                      <  2:636966/140‑1 (MQ=255)
                                                                                                                  tgCAGCATTAACGCGAATCGGATAACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTtata                      <  2:636943/140‑1 (MQ=255)
                                                                                                                  tgCAGCATTAACGCGAATCGGATAACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTtata                      <  2:636979/140‑1 (MQ=255)
                                                                                                                       cacaaaCGCGAATCGGCACACCGGGATTGAGAAACCACGCCCACAAAAACGACGCCAGATAATCACTCCGCGCACCCGCACCAGAACCAGGAAAATAATCGCCGCCCAGACAACTCATTTCTGACCGCAGTTATAACGcc                 <  2:1404298/136‑1 (MQ=255)
                                                                                                                       caTTATCGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCTGATAATCACGCCGCTCCCCAGCAGCAGAAACAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGcc                 <  2:573033/140‑1 (MQ=255)
                                                                                                                       caTTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGcc                 >  2:1140443/1‑140 (MQ=255)
                                                                                                                       caTTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGcc                 >  1:364902/1‑140 (MQ=255)
                                                                                                                       caTTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGcc                 >  1:948174/1‑140 (MQ=255)
                                                                                                                       caTTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGcc                 <  1:779467/140‑1 (MQ=255)
                                                                                                                       caTTAAAGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGcc                 <  2:571408/140‑1 (MQ=255)
                                                                                                                         ttAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCCGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCg               <  1:1147939/140‑1 (MQ=255)
                                                                                                                         ttAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCg               <  2:1357603/140‑1 (MQ=255)
                                                                                                                            aCGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCGgtg            >  2:294789/1‑140 (MQ=255)
                                                                                                                            aCGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCGgtg            <  1:604561/140‑1 (MQ=255)
                                                                                                                                 aaTCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCGGTGTCAcc       >  1:481125/1‑140 (MQ=255)
                                                                                                                                 aaTCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCGGTGTCAcc       >  1:141698/1‑140 (MQ=255)
                                                                                                                                 aaTCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCGGTGTCAcc       >  1:799102/1‑140 (MQ=255)
                                                                                                                                 aaTCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCGGTGTCAcc       >  1:1131846/1‑140 (MQ=255)
                                                                                                                                   tCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCGGTGTCACCTa     <  1:1166184/140‑1 (MQ=255)
                                                                                                                                      gATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATATCGCCCGGTGTCACCTACtt  >  2:1419184/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAGAAAGCAGGCGCAAAATAAGGACGCCAGATAATCACGCCGCTCACCAGCAGCAGAACCAGGAAAATAATCGCCGCCCAGAAAACGCATTTCTGACCGAAGTTATAACGCCCGGTGTCACCTACTT  >  NC_000913/4081289‑4081562

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