Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,707,031 C→A 6.8% P38H (CCC→CAC) ‡ rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,707,0310CA 6.8% 264.1 / 7.8 89P38H (CCC→CAC) ‡rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (35/47);  new base A (3/3);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

ACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGT  >  NC_000913/1706883‑1707179
                                                                                                                                                     |                                                                                                                                                    
acacCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTAccc                                                                                                                                                     >  2:98224/1‑150 (MQ=255)
 cacCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTAcccc                                                                                                                                                    <  1:490310/150‑1 (MQ=255)
 cacCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTAcccc                                                                                                                                                    >  1:1436308/1‑150 (MQ=255)
   ccATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTg                                                                                                                                                  >  2:218955/1‑150 (MQ=255)
   ccATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCCGCATCCATCCACCGGAGATGAAAACACAGTCCAACGGTACACCCCTCCGCCAGGTACCCCTg                                                                                                                                                  >  1:90589/1‑150 (MQ=255)
    cATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCCCCCGGTTCCCCggg                                                                                                                                                 >  1:1751042/1‑150 (MQ=255)
     aTTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTGCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGgc                                                                                                                                                <  1:527167/150‑1 (MQ=255)
        cccGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCa                                                                                                                                             >  2:636207/1‑150 (MQ=255)
        cccGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGTACCTGGCGCa                                                                                                                                             <  2:678561/150‑1 (MQ=255)
          cGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGc                                                                                                                                           <  2:416817/150‑1 (MQ=255)
           gTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGTACCTGGCGCAggg                                                                                                                                          <  1:707777/150‑3 (MQ=255)
           gTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGTACCTGGCGCAggg                                                                                                                                          <  1:833558/150‑3 (MQ=255)
            tGCGTATCATTCCCGGGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGt                                                                                                                                         <  2:192501/150‑1 (MQ=255)
             gCGTATCCTTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCGGCGATTTCTACGGCGGCAGCCATGCTCCGGAGATGAAACCCCAGTCCATCGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGtt                                                                                                                                        <  1:49948/150‑1 (MQ=255)
              cGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGttt                                                                                                                                       <  2:1692877/150‑1 (MQ=255)
                tATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTg                                                                                                                                     <  2:1486667/150‑1 (MQ=255)
                tATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTAAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTg                                                                                                                                     <  2:1393265/150‑1 (MQ=255)
                 aTCATTCCCGTGGAACACCATGCTTATGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGt                                                                                                                                    <  2:1225414/150‑1 (MQ=255)
                 aTCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCAGCGCCAGGTACCCCTGGCGCAGCGTTTTGt                                                                                                                                    >  1:969058/1‑150 (MQ=255)
                  tCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGtt                                                                                                                                   <  2:1444769/150‑1 (MQ=255)
                    aTTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTAt                                                                                                                                 <  1:1626592/150‑1 (MQ=255)
                    aTTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCATGGCGGAGCGTGTTGTTAt                                                                                                                                 >  2:159610/1‑150 (MQ=255)
                     ttCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTAtt                                                                                                                                >  2:21128/1‑150 (MQ=255)
                        ccGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTccc                                                                                                                             >  1:1217559/1‑149 (MQ=255)
                          gTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACt                                                                                                                           <  1:1194332/150‑1 (MQ=255)
                            ggAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGa                                                                                                                         >  2:1264525/1‑150 (MQ=255)
                            ggAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCACCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGa                                                                                                                         >  1:412812/1‑150 (MQ=255)
                             taaCACCATTCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACTGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGaa                                                                                                                        <  2:1453222/149‑1 (MQ=255)
                             gAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAAATAAAATCTGTGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGa                                                                                                                         <  2:874761/150‑1 (MQ=255)
                                  ccATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGc                                                                                                                   >  2:1572166/1‑150 (MQ=255)
                                   cATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCa                                                                                                                  >  1:1070747/1‑150 (MQ=255)
                                         tAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTgg                                                                                                            >  2:37619/1‑150 (MQ=255)
                                         tAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTGTTGTTATTCCACTGAAACAGCATATTgg                                                                                                            >  1:469417/1‑150 (MQ=255)
                                           aGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGgcg                                                                                                          <  1:378763/150‑1 (MQ=255)
                                            gTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCTTCCATCCACCGGAGATGAACCCCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGgcgc                                                                                                         <  2:78776/150‑1 (MQ=255)
                                              tATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTg                                                                                                       <  1:98224/150‑1 (MQ=255)
                                                  ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGTGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAgg                                                                                                   <  1:1072706/150‑1 (MQ=255)
                                                     tGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGa                                                                                                <  1:568720/150‑1 (MQ=255)
                                                      gCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAg                                                                                               <  2:1219022/150‑1 (MQ=255)
                                                      gCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAg                                                                                               <  1:869155/150‑1 (MQ=255)
                                                        aTTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTACAACGGAACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTAATTCCACAGAAACAGCATATTGGCGCTGAAGGTGAGtt                                                                                             >  1:953177/1‑150 (MQ=255)
                                                         ttCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTCGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtg                                                                                            >  1:1162926/1‑150 (MQ=255)
                                                              caaaaaaTAAACTCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCg                                                                                        <  2:705605/149‑1 (MQ=255)
                                                                 aaaataaaatCTGTGATTTCAACTGCGGCATACATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTa                                                                                     <  2:1780893/150‑1 (MQ=255)
                                                                  aaataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAg                                                                                    <  1:1074827/150‑1 (MQ=255)
                                                                      aaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGACCAACGGTACACCCCTGCGCCACGTACCCCTGGCGCACCGTGATCTTATTACACTGAAACAGCATATTGGCCCAGAAGGTCAGTTGTGCGTTAGAGta                                                                                >  1:1819509/1‑149 (MQ=255)
                                                                       aaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAACCCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCAGTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGGGTTTGCGTCg                                                                               >  1:1580672/1‑150 (MQ=255)
                                                                         atCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGc                                                                             <  1:1473093/150‑1 (MQ=255)
                                                                         atCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGc                                                                             <  1:1777782/150‑1 (MQ=255)
                                                                            tGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGAACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGAGTTAGCGTCGGCGAt                                                                          >  2:217969/1‑150 (MQ=255)
                                                                                 tttCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGt                                                                     >  1:1768595/1‑150 (MQ=255)
                                                                                    cAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTTTTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTTTGCGTTTGCGTCGGCGATAAAGTAtt                                                                  >  2:1248019/1‑150 (MQ=255)
                                                                                         gcggcATCCCTCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGcgcg                                                             <  2:90589/150‑1 (MQ=255)
                                                                                          cggcATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCGGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGAAAAAGGATTGCGCgg                                                            >  2:1568559/1‑150 (MQ=255)
                                                                                                 catccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCCGGGGCGCCGGGCGCCGCGTTTTGGTATTCCCCTGCAACAGCAGAGTGGGGGTGCAGGGGGGTTGGGGGTTGGCGGGGGGGGTCAAGGTTTGGGCGGCccgccg                                                     >  2:787736/1‑150 (MQ=255)
                                                                                                   tccaCCTGAGATGAAAACCCAGTCCAACGGTCCACCCCTGCGCCACGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTTTTGCGCGGCCAGCCGCt                                                   <  2:1579033/150‑1 (MQ=255)
                                                                                                        cGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTAccc                                              <  2:1441133/150‑1 (MQ=255)
                                                                                                         ggAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCg                                             <  1:1368008/150‑1 (MQ=255)
                                                                                                           agaTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTg                                           <  2:950754/150‑1 (MQ=255)
                                                                                                             aTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGt                                         <  2:1288466/150‑1 (MQ=255)
                                                                                                              tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTTCGTTAGCGTCGGCGATAAAGTATTTCGCGGCCAGCCGCTTACCCGTGGTc                                        >  1:855709/1‑150 (MQ=255)
                                                                                                              tGAAAACCCAGTCCAACGGTACACCCCTGCCCCCGGTACCCCTGGCGCCGCGGTTTTTTATTGCCCTGAAACAGCATAATGGCGCTGAAGGTGAGTTGTGCGTTTTCGTCCGCGATAAAGTATTTCGCGGCCACCCGCTTTCCCGTGGta                                        >  2:1370329/1‑149 (MQ=255)
                                                                                                                 aaaCCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTcgcg                                     <  1:31223/150‑1 (MQ=255)
                                                                                                                  atcccAGTCCAACGGTACACCCCTGCGCCAGGGACCCCTGTCGCAGCGTTTTGTTATTCAACTGAAACAGCATATTGGCGCTGATGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCgg                                    <  2:469417/148‑1 (MQ=255)
                                                                                                                   aCCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGc                                   <  1:583721/150‑1 (MQ=255)
                                                                                                                   aCCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTGTTGTTTTTCCACTGAAACAGCATATTGGCGCTGAAGGGGAGTTGTGCGTTAGCGTCGGCGATAAAGTCTTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGc                                   <  1:1366039/150‑1 (MQ=255)
                                                                                                                   aCCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGGCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGc                                   <  1:1358553/150‑1 (MQ=255)
                                                                                                                    cccAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCa                                  >  1:1453992/1‑150 (MQ=255)
                                                                                                                    cccAGTCCAACGGTACACCCCTGCGCCAGGTACCCCAGGCGCAGCGTTTTGTTATTCCACAGAAACAGCATATTGGCGCAGAAGGTGAGTTGTGCGTTTGCGTCGGCGATAAAGTATTTCGCGGCCAACCGCTTACCCCTGGTCCCGGCa                                  >  2:546173/1‑150 (MQ=255)
                                                                                                                      cAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCaaa                                <  1:1453195/150‑1 (MQ=255)
                                                                                                                        gTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAAt                              <  1:1284329/150‑1 (MQ=255)
                                                                                                                         tCCACCGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGGGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAAtg                             <  2:1751042/150‑1 (MQ=255)
                                                                                                                         tCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTTGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTTAGCGGCCACCCGCTTACCCGTGGTCGCCGCAAAAtg                             >  2:1314382/1‑150 (MQ=255)
                                                                                                                           cAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAAtgct                           >  2:1210753/1‑150 (MQ=255)
                                                                                                                             aCGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAAtgctgc                         <  1:1264525/150‑1 (MQ=255)
                                                                                                                              cGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGcc                        <  2:1572324/150‑1 (MQ=255)
                                                                                                                              cGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGcc                        <  1:1572166/150‑1 (MQ=255)
                                                                                                                              cGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCCTCGGCGATAAACTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGcc                        >  2:1688272/1‑150 (MQ=255)
                                                                                                                                gTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTg                      <  2:977029/150‑1 (MQ=255)
                                                                                                                                    acCCCTGCGCCAGGTCCCCCTGGCGCAGCGTGTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCa                  <  1:315481/150‑1 (MQ=255)
                                                                                                                                     ccccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAc                 >  1:61181/1‑150 (MQ=255)
                                                                                                                                     ccccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACACGTGGTCGCGGCAAAAAGCTGCCTGTTCAc                 >  2:810565/1‑150 (MQ=255)
                                                                                                                                     ccccTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGGGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAc                 <  1:1272771/150‑1 (MQ=255)
                                                                                                                                     acgcTGCGCCAGGGGTACCTGGCGCAGGGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAAAGCTGCCTGTTCAc                 >  1:1764769/4‑150 (MQ=255)
                                                                                                                                       ccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAcgc               >  2:709420/1‑150 (MQ=255)
                                                                                                                                        cTGCGCCAGGGGTACCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCCGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAcgcg              >  2:1415716/1‑150 (MQ=255)
                                                                                                                                         tGCGCCAGGGGTACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACgcgc             >  2:930995/1‑150 (MQ=255)
                                                                                                                                          gcgcCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGcc            <  1:1825898/150‑1 (MQ=255)
                                                                                                                                             ccAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCAc         >  1:697410/1‑150 (MQ=255)
                                                                                                                                                 gTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCg     >  1:288462/1‑150 (MQ=255)
                                                                                                                                                 gTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCg     >  1:813604/1‑150 (MQ=255)
                                                                                                                                                  tACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCTGCAAAATGCTGCCTGTTCACGCGCCCACCTCgg    <  1:1617541/150‑1 (MQ=255)
                                                                                                                                                    ccccTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGt  >  1:1637506/1‑150 (MQ=255)
                                                                                                                                                    ccccTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACATCGGGt  >  1:187504/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGT  >  NC_000913/1706883‑1707179

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